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The repercussions of an inaccurate account of electronic states near the Fermi level Ef by hybrid 
functionals in predicting several important metallic properties are investigated. The difficulties in¬ 
clude a vanishing or severely suppressed density of states (DOS) at Ef, signihcantly widened valence 
bandwidth, greatly enhanced electron-phonon (el-ph) deformation potentials, and an overestimate 
of magnetic moment in transition metals. The erroneously enhanced el-ph coupling calculated by 
hybrid functionals may lead to a false prediction of lattice instability. The main culprit of the 
problem comes from the simplistic treatment of the exchange functional rooted in the original Fock 
exchange energy. The use of a short-ranged Coulomb interaction alleviates some of the drawbacks 
but the fundamental issues remain unchanged. 


As one of the most successful and extensively imple¬ 
mented theories in electronic structure methods, the im¬ 
portance of density functional theory (DFT) [1, 2] can not 
be overstated. However, despite its tremendous triumph 
in describing many ground state properties, DFT within 
the local density approximation (LDA) [2, 3] or the gen¬ 
eralized gradient approximation (GGA) [4-6] has seri¬ 
ous limitations when it comes to describing excited state 
properties of materials, a problem that is often loosely 
referred to as the LDA band gap problem. 

Strictly speaking, the Kohn-Sham (KS) eigenvalues 
cannot be interpreted as the quasiparticle energies in 
solids. Nevertheless, there has been much effort put 
into the development of better energy functionals within 
DFT that can also provide a more faithful description 
of the quasiparticle energies. One class of such function¬ 
als, known as hybrid functionals [7], in which the con¬ 
ventional LDA or GGA exchange-correlation functional 
is mixed with a fraction of the Hartree-Fock (HF) ex¬ 
change energy, has the promises of providing a better 
description of a wide range of materials properties [8- 
15], including the band gap of semiconductors. There 
are notable exceptions, however. For example, the ex¬ 
perimental band gap of FeS 2 , an earth abundant mate¬ 
rial with a very high optical absorption coefficient in the 
visible region, is about 0.95eV [16, 17], but the Heyd- 
Scuseria-Ernzerhof (HSE) [18] hybrid functional predicts 
a band gap of about 2.6 eV. [19] 

We show that the application of hybrid functionals to 
metallic systems gives rise to a range of false or inaccu¬ 
rate results. These include, but not limited to, 1) en¬ 
larged valence bandwidth, 2) vanishing or greatly sup¬ 
pressed density of states (DOS) at the Fermi level, 3) 
significantly overestimated el-ph deformation potentials, 
4) false predictions of lattice instability for systems (e.g., 
MgB 2 ) with a strong el-ph coupling, and 5) enhanced 
magnetism in transition metals. The electronic structure 
near the Fermi level is one of the most important proper¬ 
ties of metals. Many fundamental response, thermal, and 


transport properties of metals are controlled or strongly 
affected by the electronic structure near the Fermi level. 
The inaccurate prediction of the strength of the el-ph 
coupling, lattice stability, and magnetism are just a few 
manifestations of the fundamental difficulties that are in¬ 
herent to hybrid functionals when they are applied to 
metallic systems. We believe that the aforementioned 
problems are serious enough to justify a thorough evalu¬ 
ation of the applicability of hybrid functionals to metals 
and degenerate semiconductors. 

Our density functional calculations are performed us¬ 
ing the Quantum Espresso [20] package. The GGA func¬ 
tional of Perdew, Burke, and Ernzerhof (PBE) [21] and 
hybrid functionals, both the PBEO [22] and the HSE [18] 
functionals, are used in this work. Electron-ion inter¬ 
actions are described by norm-conserving pseudopoten¬ 
tials [23]. The plane wave basis has an energy cutoff of 
20 Ry for the calculations of Na, 50 Ry for MgB 2 , and 
180 Ry for transition metals. For transition metals, all 
semicore subshells (i.e., 3s, 3p, and 3d) are included as 
valence electrons. 

We first investigate the electronic structure of metals 
calculated with hybrid functional PBEO using Na as an 
example. Figure 1 compares the DOS of Na calculated 
using both the PBE (blue) and the PBEO (red) function¬ 
als. The DOS are calculated using a very high fc-point 
density of 40 x 40 x 40 with a small gaussian smearing 
parameter of 0.01 eV. It is well-established that the PBE 
functional gives a satisfactory description of the elec¬ 
tronic structure for a wide range of materials including 
metals. In fact, the calculated bandwidth (3.15 eV) of Na 
agrees very well with that obtained from the free-electron 
model (shown in Fig. 1 with the black curve superim¬ 
posed on the PBE result). The PBE (and free-electron 
model) result also compares reasonably well with experi¬ 
ment [24, 25] other than a well-understood quasiparticle 
renormalization factor. [26] The PBEO functional, in con¬ 
trast, gives an overly-exaggerated occupied bandwidth of 
4.2 eV. Perhaps a more alarming issue is that the PBEO 
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FIG. 1: (color online) Calculated DOS of sodium using the 
PBE (in blue) and PBEO (red) functionals. The PBE result 
agrees well with that of the free-electron model (solid black 
curve), and the PBEO result can be well understood within 
the HE theory of free electrons (solid black curve) that shows 
vanishing DOS at the Ef (shown with the dashed line). 


functional predicts a vanishing DOS at Ep as shown in 
Fig. 1. Note that the singularity is not clearly shown in 
the figure because of a finite smearing parameter used in 
our calculation. Nevertheless, the significant dip in the 
DOS at the Ep can be clearly seen. The abnormal DOS 
calculated with the PBEO functional is a direct result of 
the HF exchange energy incorporated in the functional. 
A simple free-electron calculation using the HF theory 
(with the HF exchange-energy rescaled by a factor of 0.25 
as it is used in the PBEO functional) gives nearly identi¬ 
cal results as shown with the black curve superimposed 
on the red curve (PBEO results). 

The Fermi-surface property is perhaps one single most 
important property of metals. A wrong description of 
the Fermi-surface properties, including the DOS and the 
electron energy dispersion near the Fermi level, will have 
serious consequences when it comes to predicting other 
fundamental properties of metals. In order to demon¬ 
strate the repercussion of the incorrect description of 
the electronic structure of metals near the Fermi level, 
we investigate some electron-phonon coupling related as¬ 
pects of metals using MgB 2 , a well-established multi-gap 
phonon-mediated superconductor [27-31], as an exam¬ 
ple. Quantitative predictions of the superconductivity 
in MgB 2 , including the precise calculation of the super¬ 
conducting transition temperature Tc, is one of the cel¬ 
ebrated [29-31] successes of modern electronic structure 
methods. This success would not be possible without 
an accurate description of the electronic structure near 
Ep. Figure 2 shows the DOS of MgB 2 calculated using 
the PBEO (shown in red) and PBE (in blue) functionals. 
Similar to the case of Na, the PBEO functional gives a 
significantly widened occupied bandwidth and a greatly 
suppressed DOS near Ep. In fact, a fully-converged 



FIG. 2: (color online) Density of states of MgB 2 calculated us¬ 
ing the PBE (blue) and PBEO (red) functionals. The dashed 
line indicates the Fermi level Ep. 



FIG. 3: Schematics of the polarization vector of one of the 
i? 2 g phonon modes in MgB 2 . 

PBEO calculation (with respect to the fc-grid density) 
should give a precisely zero DOS at Ep. 

It is now well-established that the strong coupling be¬ 
tween the E 2 g phonons (shown in Fig 3) and the boron- 
derived pa states is largely responsible for the enormously 
high Tc of MgB 2 . Figure 4 shows the band structure of 
MgB 2 calculated with boron atoms displaced by 0.05 a.u. 
from their equilibrium position along one of the i? 2 g 
phonon modes shown in Fig. 3. The PBE results are 
shown with blue solid curves and the PBEO results are 
shown with red dashed curves. The degeneracy of the 
pa states along the P —A direction is lifted by the E 2 g 
phonon distortion. An and Pickett [29] used this energy 
shift (splitting) to calculate the el-ph deformation poten¬ 
tial of the i? 2 g phonons and successfully estimated the Tc 
to be in the range of 32 ^ 46 K using the PBE functional. 
Our calculated E 2 g deformation potential Dp^^ using the 
PBE functional is about 12.5 eV/A, which agrees well 
the result of An and Pickett [29]. The PBEO functional, 
on the other hand, predicts a deformation potential of 
26.3 eV/A which nearly doubles that predicted by the 
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FIG. 4: Band structures of MgB 2 calculated with a frozen- 
phonon distortion along one of the E 2 g modes using the PBE 
(blue solid line) and PBEO (red dashed line) functionals. The 
magnitude of the distortion is 0.05 a.u. 



Egg displacement (a.u) 


FIG. 5: Total energy of MgB 2 as a function of the magnitude 
of a frozen-phonon distortion along one of the E 2 g phonon 
modes calculated using the PBE (red), HSE06 (green) and 
PBEO (blue) functional. 

PBE functional. This overestimated splitting between 
the two pa bands with the presence of the E 2 g phonon 
distortion is consistent with the understanding that hy¬ 
brid functionals tend to enhance the separation between 
occupied and unoccupied electronic states. Note that the 
PBEO band structure is constructed from a self-consistent 
calculation on a 24 x 24 x 24 uniform fc-grid. 

It should be mentioned that the strength of the el-ph 
coupling in MgB 2 is so strong that the E 2 g phonons in 
MgB 2 are greatly renormalized as compared to the E 2 g 
phonons in AIB 2 . [32] In fact, the E 2 g phonons become so 
soft that the system appears to be close to a structural 
instability. The greatly exaggerated el-ph deformation 
potential calculated with hybrid functionals suggests that 
these functionals might actually predicts a structural in¬ 


stability for MgB 2 . To this end, we have calculated the 
total energy of MgB 2 as a function of the E 2 g phonon dis¬ 
placement. As shown in Fig. 5, the potential curve cal¬ 
culated using the PBE functional (shown in red) shows a 
familiar third-order anharmonicity. The potential curve 
calculated using the HSE functional (green) suggests a 
greatly softened E 2 g phonon compared with the PBE re¬ 
sult (red), and the result from the PBEO functional shows 
a double-well structure, indicating an instability against 
the £' 2 g-like distortions. 

The combination of an inaccurate description of the 
electronic structure near the Fermi level and the wrong 
predictions of el-ph coupling and lattice stability suggests 
that hybrid functionals, especially the PBEO functional, 
should probably never be used for metal systems. The 
inaccurate or false results for simple systems such as Na 
or MgB 2 predicted by hybrid functionals, however, may 
not seem to be disturbing enough since these functionals 
are usually not used for simple systems. However, one 
can also argue that a reliable energy functional should 
work equally well for both simple and more complex (e.g., 
strongly correlated) systems. It should be pointed out 
that hybrid functionals have been unbiasedly applied to 
metallic systems, and the aforementioned issues are gen¬ 
eral and will be present in other metallic systems as well. 

Finally, we examine further into the limitations and ap¬ 
plicability of hybrid functionals to the understanding of 
magnetic metals. One of the advantages of hybrid func¬ 
tionals (compared with the EDA or GGA functionals) is 
that they improve the exchange energy functional and 
reduces the self-interaction error, which helps to open 
the band gap for semiconductors. This advantage should 
be examined with great caution since it also produces 
erroneous predictions for a range of metallic properties 
as discussed above. Another consequence of including a 
fraction of the HE exchange energy (either in the origi¬ 
nal form such as that in the PBEO functional or in the 
screened form as in the case of HSE functional) is that it 
enhances the localization of electronic states and the ex¬ 
change splitting between the occupied (or majority spin) 
and unoccupied (or minority spin) states. This enhance¬ 
ment, however, does not guarantee a better description 
of magnetism. 

To illustrate this issue, we have performed spin polar¬ 
ized calculations for iron and nickel using the PBE and 
the HSE hybrid functionals. Table I compares the cal¬ 
culated magnetic moments of these two systems with ex¬ 
periment. The HSE06 functional overestimates the spin 
magnetic moment by about 29% and 61% for Fe and 
Ni, respectively. On the other hand, the magnetism pre¬ 
dicted by the PBE functional is in a reasonable agreement 
with experiment [33]. We mention that similar problem 
has also been discussed recently [34, 35] that hybrid func¬ 
tionals predict wrong magnetic ground states for Co, Ni, 
Ph, Pd and Pt. 

In conclusion, we have shown that there are serious 
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TABLE I: Calculated magnetic moments (/^B/atom) and rel¬ 
ative errors in comparison with experiment [33] for Fe and Ni 
using PBE and HSE06 functionals. 



PBE 

HSE06 

Experiment 


M 

Error 

M 

Error 

M 

Fe 

2.29 

3% 

2.87 

29% 

2.22 

Ni 

0.70 

13% 

1.00 

61% 

0.62 


limitations of hybrid functionals in predicting fundamen¬ 
tal properties of metals. The inaccurate description of 
the electron structure of metals, including a vanishing 
DOS at the Fermi level, exaggerated occupied band¬ 
width, and enhanced splitting between occupied and un¬ 
occupied states, have profound and adverse effects on the 
prediction of other important properties of metals. These 
difhculties can be traced back to the overly simplihed 
HF-like exchange energy included in hybrid functionals 
which is known to give an the inaccurate description of 
electronic structure near the Fermi level. Introducing 
a short-ranged Coulomb potential as in the case of the 
HSE functional does alleviate some of the problems, but 
major difficulties remain unchanged. Although compu¬ 
tationally demanding, hybrid functionals are becoming 
increasingly popular with the implicit assumption that 
they will give a more faithful description of the underly¬ 
ing electronic structure. Our results raise a serious ques¬ 
tion regarding the applicability of hybrid functionals to 
the study metallic systems. 
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